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Abstract :The present numerical study was performed by 
solving conservation equations of mass, momentum and 
energy along with two equation based k-e model for a 
louvered horizontal cylindrical pipe by finite volume method. 
It was found that pipe air suction increased with increase in 
opening area and number of nozzles for both pipe 
configurations. Increasing the pipe diameter made air suction 
more effective in pipe configuration B (open entrance of the 
pipe). 
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I. Introduction 

Release of hot exhaust gases (temperature of about 70 °C) 
directly into the atmosphere is causing various environmental 
problems like smog formation and constant increase in 
atmospheric temperature. If the exhaust gases can be cooled, 
then the environment can be protected from its adverse effects. 
The cooling of exhaust gases can be done by sucking fresh air 
from the atmosphere through louvers of the pipe. Some 
research work was conducted in the past both theoretically and 
experimentally to overcome the harmful effects of release of 
hot exhaust gas. Han and Mungal (2001) in their study found 
that the increase in co-flow speed reduces the entrainment for 
both reacting and non-reacting jets. They also found that the 
entrainment in the near field of reacting jets was reduced by a 
factor of 2.5 due to the heat release effect. Mishra and Dash 
(2010) in their numerical investigation found that air suction 
rate is enhanced with increase in louvers opening area but 
independent of number of louvers per row. Singh et. al. (2010) 
found in their research that the ratio of tube diameter and inlet 
jet diameter along with the ratio of free stream air density to 
density of air at jet inlet condition were the main factors that 
influenced the entrainment and mixing. Barik et al. (2014) in 
an experimental and numerical investigation on air entrainment 
in infra-red suppression device found that the entrainment rate 
was reduced by 19.2% as the ratio of funnel overlap height and 
nozzle diameter changed from 0 to 1.6. Singh et al. (2003) in an 
experimental investigation found that non circular jets provided 
greater entrainment and mixing with ambient fluid than circular 
jets. They concluded that by shifting the location of jet, the 
entrainment was enhanced by 30%. 

The work found in the literature does not site the use of CFD 
methods to determine the effect of various parameters on mass 
suction rate of air into a horizontal louvered pipe for a vehicle, 
and hence the present work is done to arrive at the optimum 
design of the pipe by investigating fluid flow and heat transfer 
characteristics. The study is carried out assuming exhaust pipe 
as the nozzle and hot air as the nozzle fluid and treated to be 



incompressible. We have attempted to match the existing 
experiment results of Singh et. al. (2003) with the present 
computational methods. 

II. Mathematical Formulation 

A cylindrical pipe of diameter D P and length L P was used to 
carry out the computational investigation. The pipe is closed at 
the entrance end and its exit end is open to the surrounding 
atmosphere for Configuration A whereas for Configuration B, 
both the ends of the pipe are open to the atmosphere. A nozzle 
of diameter D n having a protrusion length L Pt (into the pipe), is 
located at the centre of the entrance face of the pipe. A 
cylindrical computational domain of diameter D CD (about 6 
times of D P ) and length L CD (about 3 times of L P ) is placed 
around the pipe so that the boundary conditions can be applied 
on the computational domain for the air to be sucked in by the 
louvers in the pipe. The louvers are made in circular shape (D h 
= 0.03 m). A total of 12 louvers are placed in each stack with a 
distance of 0.5 m being maintained between each stack. The 
flow field in the domain is computed by using three- 
dimensional, incompressible Navier-Stokes equations and a 
standard two equation k-s turbulence model with energy 
equation. The two equation k-e turbulence model is believed to 
give the correct results for such Re flows. The fluid used for 
simulation is air at 343 K with a mass flow rate of 0.02 kg/s at 
nozzle outlet for every case and is treated as incompressible. 
Steel is declared as the solid material used for all the walls 
used for the present numerical research. 
Governing Equations 

The governing equations used for the present analysis can be 
written as: 
Continuity equation: 

dx i 

Momentum equation: 
D(pU.) dp 3 r MJ k SU, 

- l/ =-— + [M-C L + -)-pUM.\ 

Dt cbc dx, dx: 3x: 1 



P=RT 



(2) 
(3) 



The density p is taken as the function of temperature 
according to the ideal gas law as per equation (3), while 
dynamic viscosity u and thermal conductivity are kept 
constant. 
Energy equation: 

£>(pT) = d u | fi t dT (4) 
Dt dx ; Pr Pr f dx 
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Turbulence kinetic energy, k: 
-£-(pk) = A+pP-pe 
Dissipation rate of k: 



£(ps) 
Dt 



k 



-v r ( — L + — -): 



v, =0.09 



5 jU, 0<|) 

3x . ^ & . 



P = -U t Uj 



C.V: 



(5) 

(6) 

(7) 
(8) 

(9) 



o k and a g are the Prandtl numbers for k and 8. 
R is characteristic gas constant and is equal to 0.287 kJ/kg-K. 
The following constants are used in the k-s equations. 
C l8 = 1.44, C 2e = 1.92, C M =0.09, o k = 1.0, o, = 1.3, Pr t = 1. 
Although these constants are normally used for internal 
flows, Fluent's default constants have not been changed for 
the present case of mixed flows, i.e. both internal and 
external. 

Boundary Conditions 

The boundary conditions for the louvered pipe can be seen 
from Fig. 1. The pipe and nozzle walls are solid and a no-slip 
boundary condition has been imposed. Pressure outlet 
boundary condition has been given to the top, side and 
bottom surface of the computational domain. Nozzle outlet 
has been given velocity inlet and temperature boundary 
condition to supply hot air into the pipe. The velocity will be 
computed from the local pressure field to satisfy continuity at 
the pressure outlet boundary while all the other variables like 
T, k and s will be computed from zero gradient condition. 
The turbulent intensity at nozzle inlet has been set to 2%, 
while 5% has been set for back flow turbulent intensity at all 
pressure outlet boundaries. 
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Fig. 1 Schematic view of computational domain with applied 
boundary conditions 

III. Numerical Solution Procedure 

The computational domain was discretized by using 
tetrahedral cells. The three dimensional equations for mass, 
momentum and energy were discretized over control volume 
using finite volume technique to obtain a set of algebraic 
equations. These algebraic equations were solved by 
algebraic multi grid solver of Fluent 14 by iteration using the 
imposed boundary conditions stated earlier. First order 
upwind scheme was considered for discretized momentum 
and turbulent equations. The effect of domain size and 
number of grids has been discussed in detail separately in the 
results and discussions. 

The SIMPLE algorithm with Standard scheme was used for 
the pressure velocity coupling. Under relaxation factors of 
0.3 for pressure, 0.7 for momentum, 0.8 for k and 8 and 1 for 
temperature were used for the convergence of all variables. 
Fine grids were used at the nozzle entrance and louver 
openings to get better accuracy. Tetrahedral cells were used 
for the entire computational domain. Convergence criterion 
for the discretized equations for the whole field residual was 
set to 10 3 for all variables except for energy equation, whose 
residual level was set to 10" 6 . 

IV. Results and Discussions 
Comparison with other computations 

The ratio Qsuctior/Qiniet obtained from numerical computation 
using k-8 turbulence model was compared with the experiment 
from Singh et. al. (2003) and Pritchard's relation for confined 
jet as can be seen from Fig. 2. The results from the present 
CFD k-8 model match well with the analytic solution 
developed by Pritchard (1977) and experiment from Singh et. 
al. (2003). The radial profile of mean axial velocity of a wall- 
constraint jet can be seen from Fig. 3. The mean axial velocity 
u(x, y) is normalised by the centreline velocity u(x, 0) and the 
radial distance is scaled by the jet half width (yo.s). The jet half- 
width is defined as the transverse distance from the jet axis to 
the location where the mean axial velocity u(x, y) is half of the 
centreline value u(x, 0). Singh et. al. (2003) and Becker et. al. 
(1963) proposed that if the velocity profile exhibits self- 
similarity, then the axial velocity profile will be reduced into a 
single curve. From Fig. 3, it can be seen that the results from 
the present CFD k-8 model for radial profiles matches well 
with the experimental results from Singh et. al. (2003) and the 
empirical relation developed by Becker et. al. (1963). It was 
also possible to match the variation of entrainment ratio 
Qsuction/Qiniet as a function of jet location for a circular jet from 
the experimental results available from Singh et. al. (2003) 
with the present CFD k-8 model as can be seen from Fig. 4. 
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Fig. 2 Comparing entrainment ratio as a function of D P /D n for 
present CFD model with Pritchard's relation and experiment by 
Singh et. al. 
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Fig. 3 Radial profile of axial velocity of the present CFD 
model compared with empirical correlation developed by 
Becker et. al. 
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Fig. 4 Variation of entrainment ratio with jet location 
Effect of louvers opening area on air suction rate 

Fig. 7 and Fig. 8 show the effect of increasing opening area on 
air suction rate as a function of nozzle fluid temperature for 
pipe Configuration A and B. The area of each row of louver 
was kept constant at 0.0085 m 2 and keeping the area of the 
entrance face of the pipe constant at 0.03142 m 2 . For pipe 
Configuration A only louver opening area was considered 



while for pipe Configuration B area of the entrance face was 
considered in addition to louver opening area. The simulations 
were carried out by increasing the rows of louvers up to 10 
symmetrically around the periphery maintaining a distance of 
0.05 m between each row. All the other parameters were kept 
constant as can be seen in the figure. From Fig. 7 and 8, it is 
clearly seen that by increasing the opening area, the air suction 
rate increased with a decrease in pipe outlet temperature. But 
after a certain opening area (about 0.051 m 2 for Configuration 
A and about 0.09 m 2 for Configuration B) there was no 
appreciable change in air suction rate because of the fact that 
the louvers located towards the top of the pipe do not create 
enough low pressure inside the pipe for the atmospheric air to 
be sucked in. 

As atmospheric pressure was imposed on the top surface of the 
pipe, the louvers placed close to the top were not very effective 
compared to the louvers placed at the entrance half of the pipe. 
The streams of hot air from the nozzle (temperature about 70- 
80 °C) and the atmospheric air (temperature about 27 °C) got 
mixed inside the pipe and left from the top surface of the pipe. 
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Fig. 7 Effect of louver opening area on air suction rate as a 
function of pipe outlet temperature 
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Fig. 8 Effect of total opening area on air suction rate as a 
function of pipe outlet temperature 

Effect of number of nozzles on air suction rate 

Fig. 9 shows the effect of number of nozzles on air suction rate 
as the function of nozzle fluid temperature. The simulation was 



UER@2015 



Page 175 



International Journal of Engineering Research 
Volume No.4, Issue No.4, pp : 173-177 



ISSN:2319-6890)(online),2347-5013(print) 

01 April 2015 



carried out for 1, 2, 4 and 6 nozzles. When single nozzle was 
used, it was located at the centre of the entrance face of the 
pipe. A constant pitch circle diameter (D pd =0.1 m) was used 
for placing the nozzles when the number of nozzles was 
increased to 2, 4 and 6. For single nozzle, the nozzle diameter 
was taken as 0.05 m. In order to maintain a constant mass flow 
rate and total cross section area of the nozzle, the diameters for 
2, 4 and 6 nozzles were taken as 0.035 m, 0.025 m and 0.02 m. 
The simulation was carried out by keeping other parameters 
constant which can be seen in the figure. It was seen that when 
the nozzle number is increased to 6, the entrainment rate was 
increased to 1.1 times and there is sharp decrease in outlet 
temperature. As the number of nozzles increased the low 
pressure is created uniformly inside the pipe and suction rate is 
becoming more. Also there is chance of better mixing of the hot 
and cold fluid inside the pipe, so the outlet temperature of pipe 
falls for both the pipe configurations when the number of 
nozzle increased to 6. 
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Fig. 9 Effect of number of nozzles on air suction rate as a 
function of pipe outlet temperature 

Effect of change in pipe diameter on air suction rate 

Figure 10 shows the effect of change in pipe diameter 
normalised by nozzle diameter on air suction rate as a function 
of nozzle fluid temperature. The louver opening area, pipe 
length, nozzle diameter and protrusion length mass flow rate 
and inlet temperature of nozzle fluid were kept constant as can 
be seen from the figure. The pipe diameter is increased from 
0.15 m to 0.2 m. It is seen from the figure that with the increase 
in pipe diameter, there is a decrease in air suction rate for 
configuration A due to which the outlet temperature of the 
nozzle fluid increases whereas the case is opposite for 
configuration B. For pipe configuration B, the entrance face of 
the pipe is open to atmosphere, which increases the opening 
area allowing more air to be sucked into the pipe and reducing 
the pipe outlet temperature. 
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Fig. 10 Effect of change in pipe diameter normalised by nozzle 
diameter on air suction rate as a function of pipe outlet 
temperature 
Conclusion 

The conclusions that can be drawn from the present 
investigation are as follows: 

• The air suction rate increased with the increase in louver 
opening area which thereby decreased the outlet 
temperature of the nozzle fluid. From the CFD analysis, it 
became clear that maximum suction effect was observed 
when the louvers were placed at the entry half of the pipe 
and closer to the nozzle. The air suction rate strongly 
depends on the opening area for both pipe configurations 
but minimum temperature was obtained in configuration B. 

• The air suction rate is increased up to 1.1 times by using 6 
nozzles as compare to single nozzle for both pipe 
configurations. 

• For a given nozzle flow rate and opening area, it was found 
that the air suction rate decreased for configuration A 
whereas increased for configuration B when the pipe 
diameter was increased. 

Nomenclature 

A h = Louvers opening area 

A Pf = Area of entrance face of pipe 

Ci e = Constant of value 1.44 

C 28 = Constant of value 1.92 

D CD = Diameter of computational domain 

Dp = Diameter of pipe 

D n = Diameter of nozzle 

Dj = Diameter of jet 

D h = Diameter of circular louver 

D pd = Pitch circle diameter 



D k =^ 



— [(H + fiL)*, 

dxj <J k dxj 



D„=_L [(U+ *.) 



dx. 



as 

dx. 



L CD = Length of computational domain 
L P = Length of pipe 
Lpt = Protruding length of nozzle 
Ji = Jet Location 

riWet = Mass flow rate through the nozzle 
m suction = Air suction rate through louvers 
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k = Turbulent kinetic energy 
P = Pressure 
Pr = v/a 

Pr t = Turbulent Prandtl number 

Re n = Reynolds number based on nozzle, pV n D n /u 

T = Temperature 

U = Velocity 

V n = Velocity at nozzle exit 
Greek Symbols 
u = Kinematic Viscosity 
p = Density 
u = Shear viscosity 
ji t = Turbulent viscosity 
s = Dissipation rate 
a k = Turbulent Prandtl number for k 
a g = Turbulent Prandtl number for 8 
§ = Scalar variable either k or 8 
a = Thermal diffusivity 
Subscripts 
n = Nozzle 
oo = Free stream 
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